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Tins is undoubtedly the movement of tbe principal star, wbicb 
has, therefore, an annual proper motion of o //f o36 in the direc¬ 
tion of 72°. 


On the Determination of Double Star Orbits from Spectroscopic 

Observations of the Velocity in the Line of Sight. By Arthur 

A. Rambaut, M.A. 

(Communicated by Sir. R. S. Ball.) 

I have shown in a paper which I read before the Royal 
Irish Academy in May 1886 how the relative velocity of the 
components of a double star is connected with the parallax and 
the elements of the orbit. In a later paper, published in the 
Monthly Notices of the Royal Astronomical Society for March 
1890, I put the results in a more convenient form, and pointed 
out what velocities we might expect to find in the already 
known optical doubles. 1 The smallness of these velocities, 
however, seemed to show that, except perhaps in a very few 
cases, the spectroscopic observations of these stars were likely to 
prove extremely difficult. 

At the time of writing the latter paper, although Professor 
Vogel had demonstrated the duplicity of Algol , and Professor 
Pickering the duplicity of /3 Aurigce and £ Ursce Majoris , the 
results were so novel and startling that, although the idea 
occurred to me of solving the inverse problem—viz., to obtain 
the orbit from the spectroscopic observations alone—it did not 
seem to me that the observations were reliable enough to 
warrant the expenditure of the time which the solution of the 
question would require. 

The beautiful results, however, obtained by Professor 
Pickering in the case of /3 Aurigce , and published in the Fourth 
Annual Report of the Henry Draper Memorial, seem to throw a 
wholly new light on the question. 

The frontispiece of this report contains five figures, the last 
four of which are reproductions of some of his photographic 
spectra, and exhibit the nature of the data on which his results 
depend. These prove at a glance that there is no question about 
the periodic displacement of the lines, and also show that the 
amount of the displacement is capable of measurement with very 
considerable accuracy. 

It is, however, the first of these figures with which we are 
principally concerned at present. This is a curve in which the 

1 At the time of writing this I was not aware of a paper by Professor 
Niven in the Monthly Notices of the R.A.S., May 1874, dealing with the sub¬ 
ject, in which he publishes formulae almost identical with those given in the 
communications mentioned above. 
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time is taken as the abscissa and the relative velocity in the line 
of sight as ordinate. The portion of the curve given covers the 
whole of the month of December 1889, containing 27 measures 
of the displacement of the lines, and in no case does the actually 



observed displacement differ from that derived from the curve 
by more than the accidental errors of measurement. 

A glance at the curve shows that there is a periodic displace¬ 
ment occurring in about four days. To determine the period 
with greater accuracy it is necessary to take a series of observa¬ 
tions extending over a long time, and to divide the interval 
which has elapsed by the number of periods which have 
occurre.d. Dor this purpose, of course, the longer the star is 
under observation, the more accurately can the period be deter¬ 
mined. Prom the portion of the curve given in the report I find 
for the period: 

P = 3*968 days. 

We can now determine the form of the curve with greater 
accuracy. Dor if we take any arbitrary epoch as the zero of 
abscissae and subtract once, twice, . . . n times P from the 
time, we can reduce all the observations within the same period, 
and thus obtain a very large number of points through which to 
draw the curve of velocities. 

Drom Professor Pickering’s report, treating the observations 
in the way just described, I obtain the following figures :— 
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In drawing the Curve of velocities from these data there are 
two conditions which must be fulfilled, the reasons for which will 
appear when we come to determine the form of the orbit. The 
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first of these is that the area of the curve above the axis of t shall 
be equal to that below the same line ; the ^second, that the area 
contained between the curve, the axis of t and the ordinate cor¬ 
responding to maximum positive velocity, shall be equal to that 
contained between the curve, the axis of t , and the ordinate cor¬ 
responding to maximum negative velocity. 



Fig. 1. Curve of Velocities. 


With these two conditions in view, the curve shown in fig. 1 
has been drawn from the figures given above. 

I shall refer to it as the curve of velocities. 


To Determine the Form of the Orbit. 

If we take the line of sight as the axis of x , and the centre of 
the principal star as the origin of coordinates, we have, if Y is 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjoumals.org/ at Florida Atlantic University on July 21, 2015 








18 91MNRAS. .51. -316R 


March 1891. 

the relative velocity 
equation 


of Double Star Orbits etc . 319 

of the components in the line of sight, the 



and consequently 



(1) 


The following considerations will enable us to determine the 
limits between which the integral is to be taken. 

If we take the line of nodes as the axis of y\ we have x 
equal to zero when the companion star is passing through this 
line. Now I have shown in my paper, already referred to, in the 
Monthly Notices , that when the companion is in this position the 
velocity in the line of sight is greatest. Hence it follows that 
if t is the epoch corresponding to a maximum velocity, we have 

t 

Vdt .(2) 



But when the companion is passing through the line of 
nodes on the other side of the primary, we have again 0=0, and 
the velocity in the line of sight again a maximum, but in the 
opposite direction. Hence, if r' is the epoch corresponding to 
this configuration, we have 



But the expression JV^ is the area of the curve in fig. 1 

between the greatest positive and the greatest negative ordinates ; 
and since we see that this must be equal to zero, we obtain the 
second condition mentioned above for guidance in drawing our 
curve of velocities. Further, since the value of x is again zero 
at the line of nodes after a complete revolution has taken place, 
we see that we have 



o; 


or, what amounts to the same thing, the area of the curve of 
velocities, measured from the greatest negative to the greatest 
positive ordinate, is equal to zero. This, taken in connection 
with what we have already obtained, gives the first condition 
mentioned above, viz. that the area above the axis of t is equal 
to that below the same line. 

In drawing this curve and computing the area contained 
between it and the axes, perhaps the most convenient method, and 
one which is quite sufficiently accurate for our present purpose, 
considering the nature of the data, is to plot down the observed 
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velocities on millimetre paper to any convenient scale, taking 
i 00 mm to represent the whole period. We then take a sheet of 
ordinary ruled writing-paper and rule it with eleven equi¬ 
distant lines, so as to divide it up into a number of spaces, as 
shown in the table given below. For every square centimetre in 
fig. 1, therefore, we have a corresponding space in the table, 
and if we enter in each of these spaces the number of square 
millimetres inclosed by the curve in the corresponding square on 
fig. 1, we have a register of the area of the curve in a very con¬ 
venient form, which enables us to see where the shape of the 
curve is in error, and to correct it as may appear necessary, so as 
to fulfil the two conditions mentioned above, without the trouble 
of going back and recomputing the whole area of the curve at 
each alteration. 



Hitherto we have only been concerned with the abscissae, but 
in addition to the equation (1) we have the relation 
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dy dx , , 

x -f- —y— = c, a constant, 
at, at . 

which may be written in the form 

*- # *I b+ J£] •••••••• .(3> 

This equation will enable ns to determine y if we have any 
means of obtaining the values of the constants c and B. 

In order to determine B it is necessary to refer to the equa¬ 
tion 

' la sin y' 

* v= ^ cos A + cos ( e ~ .(4) 


which I have obtained in my previous paper. 

As I have pointed out there, if 6 — A=o or 7r, i.e. if the body 
is passing through the line of nodes, the velocity in the line of 
sight is a maximum. Let these two values of Y be Y x and Y 2 . 
We thus obtain 


|(Y+v 2 ) 


la sin y 
P 


e cos A.. 


But when ?/=o the line joining the two components— i.e. the 
radius vector—must be at right angles to the line of nodes, and 
consequently cos ( 6 — A)±=o. Hence, if Y 3 is the velocity in the 
line of sight when in this position we have 

v 3 =J(Y+v 2 v 

or, the ordinate y is zero at the moment when the velocity is a mean 
between the greatest positive and the greatest negative velocities . 
Denoting this epoch, which we can determine from the curve of 
velocities, by the letter t q we have 

y ^f d i .(s) 

T 0 

We have not the means of determining from the spectro¬ 
scopic observations what is the Value of c, but if we take it to 
be any convenient factor, we shall obtain proportional values of y 
from equation (5), which, in connection with the values of x 
derived from equation (2), will enable us to construct an ellipse 
which is an orthographic projection of the real orbit on a plane 
passing through the line of nodes, although the angle between it 
and the real orbit is indeterminate. 

How this being an orthographic projection, the centre of the 
projected, is the projection of the centre of the real, ellipse, and 
the ratios of the various parts of the major axis of the real orbit 
are unaltered by projection. Consequently, if C (fig. 3) is the 
centre of the ellipse obtained as described above, P the primary, 
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which is also the origin of coordinates, and A and B the points 
in which the line CP meets the ellipse, we have 

CP 


CA 


— e 


( 6 ) 


We have also from equation (4) the relations 

la sin 7 v , „ ?asin7 , 

Vl = (ecos A+ i), and V 2 = — (« cog A -1); 


P-tta/i — e z 


whence we obtain 

= .(7) 

from which X may be determined. 

Again, if a represents the semi-major axis in miles, and B 
the radius of the earth’s orbit, we have 


a _a 

R tt’ 

and equation (4) becomes 

PRVvT=^ 
l (e cos A + cos (0 — A) ) 


( 8 ) 


from which a sin y may be determined, all the quantities on the 
left-hand side of the equation being already known. It will, 
however, be sufficient in general to take V l and V 2 from which 
to determine a sin y. Thus we have 

PRVi Vi — e 2 = la sin y (e cos A + 1) ; 

PRV 2 Vi — e 2 - Id sin y (e cos A — 1); 

whence 

a sm 7 = --—-(Vj-V 2 ).(9) 


There still remains one element which the data at our 
disposal enable us to determine. This is T, the epoch of peri- 
astron passage. We must, however, first determine the epoch of 
the passage through the node, which we shall call c, and which 
may at once be read off the curve of velocities in fig. 1, being, in 
fact, the time of maximum velocity. Knowing the eccentricity 
we can then find the mean anomaly nt corresponding to the 
true anomaly X. We have also got the mean motion equal to 

2 - , whence we obtain 
P ’ 

T = e±-*-£..(IC) 

2TT 


In the case before us, from the velocities given above and 
the curve in fig. 1, which is drawn by means of them, we have 
the time of greatest velocity (i.e. e) equal to 0 02 xP later than 
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the epoch which we have taken as the zero of abscissae. The 
latter epoch is December 2*65, and since o*o2xP=o*o8 days, 
we have 

e = Dec. 273. 

Again, measuring the area of the curve from the ordinate 
corresponding to this epoch, we have the following values of 

1000 

x, -, and 


t 

X 

1000 

X 

/ 1000 

\ X 

o-o 6 

+ 330 

+ 3*030 

+ 918 

•08 
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•i6 
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1*18 

*i8 
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1*021 

1*04 

*20 
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0-97 

*22 

1032 

0*969 

0*94 

•24 

1030 

0*971 

0*94 

•26 
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0*992 

0*98 

•28 

966 

i*035 

1-07 

•30 

905 

1-105 

1*22 

•32 

828 

1-208 

1-46 

*34 

735 

1-361 

1-85 

•36 

629 

1-590 

253 

•38 


1*957 

■3-83 

•40 

+ 385 

+ 2-597 

+674 

0*50 

- 310 

- 3-226 

+10-41 

*52 

448 

2-232 

498 

*54 

582 

1-718 

2-95 




t 

X 

1000 

X 

(~y 

0*56 

- 711 

— 1-406 

+1-98 

•58 

833 

1-200 

1*44 

*60 

945 

1*058 

1*12 

•62 

1045 

0*9*5 7 

0-92 

•64 

1131 

0*884 

0*78 

•66 

1204 

0831 
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•68 

1264 

0*791 

063 

•70 

1310 

0-763 

0-58 

72 

1340 

0-746 

0'56 

*74 

1354 

0739 

o*55 

•76 

1352 

0740 

0*55 

•78 

1334 

0-750 
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•80 
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0771 

0-59 

•82 

1243 

0-805 

0-65 
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*86 
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•88 
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■90 
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i*i44 

I-3I 

•92 
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i*337 

1-79 

*94 
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+ 2*71 


It will be noticed that the values of x from £=0*40 to 
t= 0*50 and from ^=0*94 to t=o‘o 6 have been omitted from this 
table. This has been done in consequence of the fact that, x 

becoming very small, the corresponding values of .- IOO ° and 

x 

( IOQ ° ) become too large to be manageable; in fact the curve 

in fig. 2 has a pair of infinite branches corresponding to these 
epochs. 
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The next step in the process is to plot down as ordinates the 

values of ( IQQ -) which are given in the last column of this 
' x ' 

table, taking the corresponding values of t as abscissae. We 
thus obtain the interpolating curve for 



shown in fig. 2. To find the epoch r 0 it is necessary to revert to 
fig. 1, and we find that Y^ + 84 and Y 2 =— 71, whence 



Y 3 =+ 6 * 5 , and the epochs corresponding to this value of Y are 
/=o # 2i5 and t= 0*761. These two epochs are marked in fig. 2 
with double-headed arrows. Calculating the areas contained 
between these ordinates, and the ordinates corresponding to the 
epochs o’02, o'04, &c. we get the following values for 
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j 




r t 

t 

c dt 

t 

c dt 
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Since the value of c is 

perfectly arbitrary, 

we may 

for con- 


venience divide these numbers by ioo. 

We thus obtain, by means of equation (5), the following 
corresponding values of x and y :— 


t 

X 

y 

t 

X 

y 

o‘c6 

+ 3 * 3 ° 

- 5-38 

*56 

- 7*11 

+ 5*62 

•08 

4*85 

5 *i 4 

*58 

8*33 

5* 2 5 

•10 

6-24 

471 

•60 

9*45 

4*82 

12 

7‘45 

4-17 

•62 

10*45 

4*28 

*14 

845 

3 ; 46 

•64 

11*31 

3*73 

16 

9‘22 

2*67 

•66 

12*04 

312 

•18 

979 

176 

•68 

12*64 

2*40 
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3 b 
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With these values of x and y the positions marked in fig. 3 
have been plotted down, and the “apparent” ellipse drawn 



Fig. 3. “ Apparent ” Ellipse. 


throngh them. Finding C the centre of this ellipse, and draw¬ 
ing CP to meet the curve in the point A, we have 


CP 

Cl = 01 5^ = « 


by equation (6). 

We have also Y 1 +Y 2 =i3 and Yj — Y 2 = i55 : 
equation (7) 


e cos X 


J 3 

155 


0-0839. 


and substituting the value just found for e, we find 


* = ± 57 °' 43 - 


hence by 


In computing a sin y, it must be remembered that the unit 
in terms of which we have expressed Y hitherto is perfectly 
arbitrary. But when the displacements of the lines are expressed 
in wave-lengths we can at once convert the values of Y into miles. 
In the Henry Draper Report the wave-lengths are not given, but 
Professor Pickering mentions that the greatest velocity amounts 
to 150 miles per second. This would make \ (V 1 — Y 2 )==i39’5 
miles per second. Hence we obtain, by means of equation (9), 
—remembering that in that equation P is expressed in years, and 

that consequently we must put for its value— 

a sin y — 0*081 x R. 


If we assume It to be 93,000,000 miles, this corresponds to a 
distance of 7^ millions. What the actual value of a may be 
the spectroscope does not enable us to determine, since in most 
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cases the value of y will remain wholly indeterminate. In 
some special cases, however, such as that of Algol , where the 
variation of the lines synchronises with a variation in the 
brightness of the star, we may assume that y is very nearly 90°, 
and that accordingly a sin y may be taken as the actual value of 
the semi-axis major. 

We have still to compute the value of T. For this purpose 
we have the equation 

tan 1 u - A / tan £ 6 , 

2 V i + e 

whence we obtain u=5°°'2, when 0 =\. Also the mean anomaly 
corresponding to this is nt =.Hence we obtain, by means 
of equation (10) 

T = Dec. 2 65 + 0 02 . P + o*i2 . P. 

= Dec. 2*65 + 0*14 . P = Dec. 3*21. 

We thus have the following elements of the orbit of 
/3 Aurigce :— 

P = 3*968 days; e — 0*156; 

T = 1889, Dec. 3*21; \ = 57°'43 ; 

and 

a sin y = 0*081 x R = 7,500,000 miles ±. 

In order to test how far these elements will give results 
agreeing with the data from which they are derived, I have 
computed the velocity at intervals of a tenth of a period from 
the epoch of periastron passage. This occurs at /=o*i4, as we 
have already seen. 

I accordingly obtain the following table of results :— 


t 

Vc 

v 0 

t 

V c 

v 0 

0-04 

+ 83*0 

+ 83-0 

0-54 

— 60’2 

-657 

'14 

+ 48 2 

+ 45*2 

•64 

-352 

— 40-0 

•24 

-142 

- 7*0 

'74 

~ 1 '3 

- 2-3 

'34 

- 58*5 

-510 

•84 

+ 377 

+ 38*0 

*44 

-71*0 

- 7 o *5 

•94 

+ 72*1 

+ 72*2 


in which V c denotes the computed velocity, and V 0 that read 
from the curve in fig. 1. The points corresponding to the values 
of Y c are marked in fig. 1 by small crosses, from which it will be 
seen that on the whole fc they agree fairly well with the curve of 
velocities. 

It will also be observed that the greatest discrepancies occur 
in what I have called the negative part of the curve. [I may, 

3 b 2 
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however, remark in passing that positive and negative velocities 
are only relative terms in cases like the present, where the lines 
in the spectrum are of sensibly equal brightness, since the sign of 
the velocity will wholly depend on the choice we make between 
the lines as to which shall be considered as belonging to the 
primary and which to the satellite.] That these discrepancies 
should occur is hardly to be wondered at considering the scarcity 
of observations to guide us in drawing the curve in this region. 
In the positive part, where the observations are more numerous, 
the shape of the curve seems to have been hit off with very con¬ 
siderable accuracy, as the computed velocities lie very close to it 
indeed. A glance, also, will show that a very small alteration in 
the form of the negative portion of the curve would bring it into 
conformity with the computed points without doing much 
violence to any of the observed data, and I have no doubt that a 
curve might easily be drawn which would agree as well with 
the observed velocities as that shown in fig. 1, and at the same 
time give more consistent results. As, however, the velocities 
which I have used in constructing my curve are obtained by 
measuring the ordinates in Professor Pickering’s report, and 
are, in consequence of the small scale on which his curve is 
given, affected with considerable errors, it does not appear 
advisable to push the calculation any further till the numerical 
results obtained by Professor Pickering are available. I have 
lately heard from Professor Pickering to the effect that the defini¬ 
tive measures upon the photographs of £ Ur see Majoris and 
ft Aurigoe are in progress, and that he hopes before long to 
print their results with conclusions derived from them with 
regard to the orbits of the components of these stars. When 
these results appear I hope to be able to obtain more accurate 
values for the quantities P, T, e , A, and a sin y by the method 
pointed out above. In the meanwhile, perhaps, the results 
given in this preliminary notice, and the method of treating the 
question, may not be without value to those interested in the 
subject. 

One point which my examination of the question seems to 
indicate very clearly is that the observations should be taken 
not so much at equal intervals of time, but as much as possible 
so that the points derived from them may lie equally thick along 
all parts of the curve of velocities, or, in other words, so that all 
parts of the curve may be determined with equal accuracy. In 
fig. 1 it will be seen that the majority of the observations are 
grouped about the positions of positive and negative maxima, 
whereas my results show that it is of equal importance to fix the 
exact time of the coincidence of the lines as well as the epoch at 
which the velocity is an arithmetic mean between its two extreme 
values. 

I may remark, in conclusion, that in determining the orbit of 
stars of the Algol type, in which one of the components only is 
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bright enough to form a spectral image on the photographic 
plate, and in which consequently its displacement has to be 
determined by comparing it with an artificial line, the treatment 
of the observations will be somewhat different. In this respect, 
stars such as Sirius , for instance, will come under the same 
head, if we suppose the variations in their velocity in the line 
of sight to be due to orbital motion. In this class of star the 
whole velocity in the line of sight is measured, as Professor 
Yogel has done so skilfully in the case of Algol , and this is com¬ 
posed of two parts, viz: the orbital motion relative to the 
centre of gravity of the system, and the absolute motion of this 
point, the latter of which may be assumed to be constant in 
magnitude and direction. 

If we have sufficient and accurate enough observations to 
determine the curve of velocities with certainty, it will be 
possible to separate the orbital motion of the brighter component 
from the absolute motion of the system by the first of the two 
conditions mentioned above on p. 318. For if we plot down the 
velocities as in fig. 1, the curve may or may not cut the axis of 
t , which represents the zero of total velocity, according as the 
orbital is greater or less than the absolute motion. But 
we have seen above that the line representing the zero 
of orbital velocity must cut the curve, so that the areas above 
and below it may be equal. Accordingly, if we draw a line 
parallel to the axis of t. to satisfy this condition, the distance of 
this line from the axis of t will represent the absolute motion 
of the system, and the distance from it to any point on the curve 
will represent the velocity due to the orbital motion in the corre¬ 
sponding position. 

Although the materials are almost ludicrously inadequate to 
such an investigation, I have attempted to obtain some results 
from the six spectrographic observations of Algol published by 
Professor Yogel in the Astronomische Nachrichten , No. 2,947. I 
find that the motion of the centre of gravity is about —0*7 or 
—o*8 Grerman geographical miles—Professor Yogel finds —0*5 
on the hypothesis of a circular orbit—and that the eccentricity 
is something between \ and ^. 

It must, however, be understood that these figures are but 
little better than mere guesses, on account of the paucity of 
observations, notwithstanding the great accuracy and value of 
those which are available. 

January 7, 1891. 


Additional Note . 

When the lines are equal in brightness, as in the case of 
fi Awrigoe , there is an ambiguity with regard to the sign of A, and 
consequently, at first sight, with regard to the sign in equation (10), 
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since there are no means of determining whether the velocity is 
positive or negative in any given position: it is, in fact, positive for 
one component and negative for another, and there is no means of 
discriminating between the two. If, however, one of them is fainter 
than the other, we may confine our attention to the former, and 
we are then able to state that if the displacement is towards the 
red end of the spectrum the velocity is positive ( i.e . the body is 
receding from us), and conversely. 

In any case, however, the expression given for e cos A, given 
in equation (7), is positive if Yj is greater than V 2 . Conse¬ 
quently, if we take e to be the time of passage through the node 
at which the velocity is greatest (which we may call ^), A will 
lie in the first or fourth quadrant, but there is no means of 
deciding, in the case of lines of equal brightness, which of these 
values for X we are to choose. A positive value for X with a 
rotation in one direction will give the same results as a negative 
value for X with a rotation in the other direction. The question 
is, therefore, whether the motion is from periastron to or 
from & to periastron. The following considerations will enable 
us to decide this point:— 

The times of maximum velocity and the two epochs (r G ), at 
which the velocity is a mean between its extreme values, mark 
four points on the orbit, the radii-vectores to which divide the 
orbit into four sectors, each of which contains a right angle. It is 
easy to prove that the periastron must fall in that sector whose 
area is least. We thus see that the time of periastron passage will 
fall in the shortest of the four intervals into which the period is 
divided by the four epochs, marked by the two maximum velo¬ 
cities and by the velocity which is a mean between its two 
maximum values. This shortest interval will always be either 
that immediately following or immediately preceding e. In the 
former case we take the upper, in the latter we take the lower, 
sign in equation (10). 
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Double Star Measures at Windsor, New South Wales, 1889 and 
1890. By John Tebbutt. 


This communication contains the results of my double star 
measures during the years 1889 and 1890. The Grubb 8-inch 
equatoreal was employed throughout. The column headed 
“ Hour-angles ” gives the hour-angles between which the 
measures were made, and the last column the weights assigned 
from a consideration of the conditions under which the observa¬ 
tions were made: 1 denotes an unusually bad condition and .5 
an unusually good one. In the observations numbered 1, 6, 7, 
24, 25 the components were estimated to be exactly equal. 
Observations 20 and 23 were made before sunset, and 2, 3, 4, 8, 
13, 14, 21, 22, 29 in twilight. In observation 12 the driving- 
clock did not act well, and the measures were very inconsistent. 
The rapid angular motion of Lalande 4219, as indicated by the 
Washington, Leander McCormick, and Pinner observations, is 
not confirmed by the measures obtained here. The orbit of 
y Comice Australis is at present far from being satisfactorily 
determined. Mr. Po well, in his paper published in the Monthly 
Notices for March last, states that the Windsor measures of 
distance are considerably greater than they should be. A prac¬ 
tised observer will, however, be satisfied, by a mere glance at 
the components, that the distance is much greater than that 
shown by Mr. Powell’s orbit. Mr. Burnham, in Astron. Nachr. 
Band 124, p. 74, gives the distance for 1889*41 as i // *79. The 
following errata occur in my Double Star Results in the Monthly 
Notices for November 1889 :—Page 25, for Lac. 5147 read a' Crucis 
and Lac. 5147 ; page 26, for 1888*320 read 1888*326. 
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